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Abstract 

Using the contraction of the SU(3) algebra to the algebra of the rigid rotator in the large 
boson number limit of the Interacting Boson Approximation (IBA) model, a line is found 
inside the symmetry triangle of the IBA, along which the SU(3) symmetry is preserved. The 
(j [ line extends from the SU(3) vertex to near the critical line of the first order shape/phase 
transition separating the spherical and prolate deformed phases, and lies within the Alhassid- 
—i' Whelan arc of regularity, the unique valley of regularity connecting the SU(3) and U(5) vertices 
■ amidst chaotic regions. In addition to providing an explanation for the existence of the 
arc of regularity, the present line represents the first example of an analytically determined 
approximate symmetry in the interior of the symmetry triangle of the IBA. The method is 
applicable to algebraic models possessing subalgebras amenable to contraction. This condition 
is equivalent to algebras in which the equilibrium ground state (and its rotational band) 
become energetically isolated from intrinsic excitations, as typified by deformed solutions to 
. the IBA for large numbers of valence nucleons. 

PACS numbers: 21.60.Fw, 21.60.Ev, 21.10.Re 

X: 

o5 : 1 Introduction 

The study of chaotic properties of the Interacting Boson Approximation (IBA) model pQ, both 
classically and quantum mechanically, led to the discovery [21 [3] of a narrow strip of nearly 
regular behavior inside the symmetry triangle [I] of the IBA, connecting the U(5) and SU(3) 
limiting symmetries, in addition to the regular region along the U(5)-0(6) leg of the triangle. 
While the existence of the latter is known to be due to the underlying SO (5) symmetry, a 
common subalgebra of both U(5) and 0(6) present throughout the U(5)-0(6) leg [5], the 
origin of the former, called the Alhassid-Whelan (AW) arc of regularity, has remained an 
open question. 

The existence of a nearly regular region connecting U(5) and SU(3), has been corroborated 
by studies of the wave-function entropy in the IBA (6] . From the empirical point of view, it 
has been found [7j that the line corresponding to the degeneracy of the f3\ and 71 bandheads 
(2+ = 0^) closely follows the arc of regularity, and 12 nuclei closely exhibiting this behavior 
have been located. Recently, it has been realized [8] that the locus of the 2+ = 0^ degeneracy 
closely follows the line of change of stability of 7- vibrations at low energies. 
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The presence of (near) regularity presupposes the existence of some underlying (approx- 
imate) symmetry. The degeneracy of the Pi and 71 bands is a well known hallmark of the 
SU(3) symmetry of the IBA [lj. It has recently been found [9] that imposing the 2+ = 2t 
degeneracy in the IBA framework leads to a line inside the symmetry triangle of the IBA 
which, in the region between the SU(3) vertex and the shape/phase coexistence region pU] 
(separating spherical from prolate deformed shapes), is located very close to the arc of regu- 
larity, while at the same time the degeneracies predicted by SU(3), not only for the Pi and 
71 bands, but also for bands belonging to higher irreducible representations (irreps) of SU(3), 
are preserved to a very good extent [9]. This result extends the notion of quasidynamical 
symmetry (QDS), originally introduced [HJ [El EE1 Q31 EE] for describing the persistence of 
limiting symmetries along the U(5)-0(6) and U(5)-SU(3) legs of the IBA, to the interior of 
the triangle. The domain of validity of the SU(3) QDS inside the IBA symmetry triangle 
has also been considered by mean field techniques [16]. The analysis of Ref. [9] was limited 
to the low-lying part of the spectrum, while regularity amidst chaoticity has been discovered 
by Alhassid and Whelan through the study of the whole spectrum [21 E] • The properties of 
high-lying rotational bands built on axially deformed ground states have been studied recently 
[T7] in the IBA, showing signatures of a SU(3) QDS extending to the highest part of the IBA 
spectrum. 

In the present work, we study Hamiltonians that, in the large boson number limit of the 
IBA, approximately commute with the SU(3) generators. In addition, we derive an analytic 
expression for the locus of these Hamiltonians inside the IBA symmetry triangle. Moreover, 
this locus in fact corresponds closely to the arc of regularity. 

This represents the first analytical identification of an approximate symmetry within the 
symmetry triangle of the IBA. The proof takes advantage of the well known contraction [18] 
of the SU(3) algebra to the [R 5 ]SO(3) algebra [T91 [20] of the rigid rotator [21]. Furthermore, 
using the contraction of 0(6) to the [R 5 ]SO(5) algebra [221 [23] of the 7-unstable rotator, we 
prove that no line related to the 0(6) symmetry exists within the triangle. 

The IBA Hamiltonian used is described in Section 2, while in Section 3 the equation of the 
line corresponding to the SU(3) symmetry is derived. The SU(3), 0(6), and 0(5) symmetries 
are considered in Sections 4-6, while an alternative parametrization is presented in Section 
7. The conclusions and outlook are presented in Section 8. The commutation relations and 
matrix elements needed for the derivations are given explicitly in Appendices 1 and 3, while the 
details of the SU(3)— )-[R 5 ]SO(3) and 0(6)— )-[R 5 ]SO(5) contractions are given in Appendices 2 
and 4 respectively. 

2 The IBA Hamiltonian and symmetry triangle 

The IBA Hamiltonian used by Alhassid and Whelan [21 [3] reads 



H (r), x) = Hi + H 2 = c r]h d + 



V 



- 1 A 




(1) 



N 
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where N is the number of valence bosons, c is a scaling factor, H\ and H2 denote the first and 
the second term respectively, 

h d = d ] -d = V5(dd) {0) , (2) 

Qxl = tfd + ds)® + X {d ] d)f. (3) 

The above Hamiltonian contains two parameters, r] and \i with the parameter rj ranging from 
to 1, and the parameter \ ranging from to —y/7/2 = —1.323. The U(5) symmetry limit 
corresponds to 77 = 1, the SU(3) limit to 77 = 0, x — ~V7 /2, and the 0(6) limit to 77 — 0, % = 0. 
These symmetries are placed at the vertices of the symmetry triangle [3] of the IBA, shown in 
Fig. 1(a). In the symmetry triangle, the narrow coexistence region [10J surrounding the critical 
line [21] separating the spherical phase from the prolate deformed phase, corresponding to a 
first-order shape/phase transition [25], is also shown. It corresponds to 77 ~ 0.8. The point at 
which the critical line reaches the U(5)-0(6) side of the triangle is known to correspond to a 
second-order shape/phase transition |25j . 

Note that an alternate parametrization of the Hamiltonian of Eq. (CTJ is often used. We 
will discuss the present results in the context of that parametrization in Section 7. 



3 The SU(3) symmetry 
3.1 Commutation relations 

The SU(3) p] algebra is generated by the angular momentum operators 

U = vioMJ, (4) 

and the quadrupole operators 

Qsu(3U = tfd + rfs)? - ^d)f. (5) 

In order to have an underlying SU(3) symmetry, the Hamiltonian of Eq. ([1]) has to commute 
with the generators of SU(3). It does commute with the angular momentum operators by 
construction, since it is a scalar quantity. We will examine the special conditions under 
which the Hamiltonian also commutes (approximately) with the quadrupole operators. The 
commutation relations needed for this task are listed in Appendix 1. 

The first term of the Hamiltonian gives 

[Hi, Qf um , v ] = cv[n d , gl (3) ,j = ^{dU - s^d)f\ (6) 

Using 

in the second term of the Hamiltonian one gets the intermediate result 
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+ { X + lP) { [(rft ^ 2) ' ^u^MU + Qxh dj d)% QS(3),J}} • (8) 

In order to obtain the conditions for which the Hamiltonian of Eq. ([I]) commutes with the 
generators of SU(3), we exploit a simplification of Eq. flB]) that occurs in the large N limit. 
In this limit, the eigenvalue expression for the second order Casimir of SU(3) [see Eq. (T4TJ1) ] 
reduces to just the A 2 term for SU(3) irreducible representations (irreps) (A,/i) with A >> /i 
and hence the ground state band [which belongs to the (2N,0) irrep] becomes energetically 
isolated from all other excitations. That is, SU(3) effectively reduces to a simple rigid rotator. 
This situation is formally known as the contraction of SU(3) to R 5 [SO(3)] [19j [20] and occurs 
when the Qsu^ operators can be replaced by mutually commuting quantities [For a detailed 

explanation, see in Appendix 2 the discussion leading to Eq. Q52J)]. If the Q^jr 3 \ operators can 
be approximated by mutually commuting quantities, Eq. (JS]) greatly simplifies. These results 
are discussed in detail in Appendix 2. 

In the large N limit, where contraction occurs, the commutators in the first two terms in 
Eq. (jHJ) will vanish. Furthermore, the vanishing of the commutator 



Qsu(3)J = -r(2*2i/|l£ + vWffiL (9) 



in this limit, implies that terms containing (cftdp 1 ' can be ignored. This fact can be understood 
qualitatively as a consequence of the relevant dominance of s bosons over d bosons within the 
ground state band, especially for relatively low-lying states in the large boson number limit. 
Eq. (jSJ) can be rewritten, without using any approximations yet, as 



+ [x 



(2h(2)l 




2 2}[«^ (fe) ^ 2) )S 2) -(Q? ) (^) (fc) )l 2) ]- (10) 

In the large N limit the terms containing (d^dp k > (in the first line and in the third line) can be 
omitted. Furthermore, in the second line, can be replaced by Q*su(3)> s ^ nce ^ as seen from 
Eq. (Cj), they differ by terms (d^d)^ 2 \ which are small. In addition, in this limit Qgu(3) can 
be replaced by the intrinsic quadrupole moment (a scalar), which is Ny/~2 in the present case 
(see Appendix A2). This replacement will be justified in detail in subsection 3.2 . This result 
is perhaps familiar in the context of the well-known property of SU(3) that B(E2 : 2^ — > Of) 
goes as iV 2 [I], that is, the collectivity of yrast transition strengths increases quadratically 
with boson number. Then in the large iV limit one is left with 

[Q? ] ■ Q?\ <?S(3) J = 2v/2iV ( x + (dts - Ji)V>. (11) 
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Then in the large N limit the commutator for the second part of the Hamiltonian reads 



[H2, QS (3)) J = c(r) - 1)2V2 (x + ^ (dKs - s^d) 



(2) 
v • 



(12) 



In order to get a vanishing commutator, the coefficients of [Ss — s^d)^ 1 in Eqs. ([6]) and ( TT2|) 
should cancel, leading in the large N limit to the condition 

X(V) = A=T JL ^-—- (13) 
AW; 2^2 (1 -77) 2 V ; 

When x is taking values between — yf/2 and 0, the parameter rj takes values between 1 and 
0.789 . From the formulae reported in Refs. [7] and [25] it is clear that the critical line in the 
large N limit corresponds to r\ crit = 0.8 for x = and to r\ crit = 9/11 = 0.818 for x — ~ V^/2. 
Thus the line described by Eq. f|T3|) cannot reach the critical line, confined in the region 
between the critical line and the SU(3) vertex. 

It should be noticed that the arc of regularity found in Refs. [21 [3] has been approximately 
described by [6] 

X(v) = 2 V ~ y. (14) 

The similarity between the lines described by Eqs. (fT3~j) and ( Il4|) can be seen in Fig. 2. 
Indeed, the two equations give very similar predictions for values of 77 between and 0.6, i.e., 
from the SU(3) vertex until quite close to the critical line. 

The symmetry triangle of IBA in the Alhassid-Whelan parametrization is shown in Fig. 1(a), 
together with the arc corresponding to Eq. (fl4"l) and the line of Eq. (fT3|) . The degeneracy line 
corresponding to E(2t) = E(2+), found in Ref. is also shown (on the right of the critical 
line) for comparison. 

We see that the present line remains very close to both the E(2p) = E(2+) degeneracy line 
and the original arc line from the SU(3) vertex until quite close to the critical region, where 
both the E(2g) = E(2+) degeneracy line and the present line turn upwards, avoiding to meet 
the critical line. 

In Figs. 1(b) and 1(c), the v and A diagrams are reproduced from Ref. [3J, with the 
lines of Fig. 1(a) plotted on them. We see that the present line remains within the valley 
corresponding to the arc of regularity for most of the way from the SU(3) vertex towards the 
critical line, turning upwards a little before reaching the critical line. 

It should be noted that the present study is greatly facilitated by the fact that the position 
of the arc of regularity appears to be practically independent of the number of bosons, as 
already remarked in Refs. [21 [3J and corroborated in Ref. [5] ■ 

In summary, we have achieved by now two goals. 

1) To prove analytically the existence of a line in the parameter space of the IBA, along 
which the Hamiltonian approximately commutes with the SU(3) generators in the large N 
limit. 

2) To prove that this line closely follows the Alhassid-Whelan arc of regularity in the region 
between the SU(3) vertex and the critical line of first order shape/phase transition. 
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3.2 Matrix elements 

In going from Eq. (TIP]) to Eq. ( TTTj) . we replaced the quadrupole operator Q S ij^ by the intrinsic 

quadrupole moment. To justify this, we consider here the matrix elements of [Hi, Q S u(3) J 

and [H2, Qsu(3) J within the ground state band. We will examine the conditions under which 
these matrix elements lead to a vanishing result. In this calculation the intrinsic quadrupole 
moment appears naturally when calculating the matrix elements of the commutators of the 
relevant parts of the Hamiltonian with the quadrupole operator, and not as a result of an 
approximation, as in the previous subsection. 

Using the standard formalism for treating matrix elements of the tensor product of two 
tensor operators and the single-boson matrix elements given in Table 1, one finds for the first 
term of the Hamiltonian 

( [N] , (2N, 0) , x = 0, L 1 1 [H x , (3) ] 1 1 [N] , (2N, 0) , * = 0, L) 

2 



cqNR 1 V2L + T, (15) 



3V7 
where 

_ (2N - L)(2N + L + 1) (2L - 3)(2L + 5) 
1 (2JV-1)(2JV) (2L- l)(2L + 3)' { } 

and x is the Vergados quantum number [27], not to be confused with the parameter x of the 
Hamiltonian ([T]). It is worth remarking that the terms d's and s*d give equal contributions of 
the same sign to the final result, despite the fact that they appear in Eq. (jBj) with opposite 
signs. 

Following the procedure described in Appendix 3, one finds for the second term of the 
Hamiltonian, including in the large iV limit only the terms appearing in the second line of Eq. 
^ 2) byQi (3) 



(fTOl) after replacing by Q{ 2) 



( [N] , (2N, 0) , x = 0, L\ \ [H 2 , Q§ (3) ] 1 1 [N] , (2N, 0) , * = 0, L) 



1 J V7 



7y _c(l -V)[X + \j q R2V2LTl, (17) 
where 

"4 (2N - L)(2N + L + 1) (2L - 3) 2 (2L + 5) 2 
3 (2N -l)(2N) (2L- 1) 2 (2L + 3) 2 



(2N-L- 2)^{2N - L)(2N + L + 3) ^(2L - 1)(2L)(2L + 2)(2L + 4)(2L + 6)(2L + 7) 
(2N-1)(2N) (2L + l) 2 (2L + 3) 2 



(2A^ + L - 1)^(2^ - L - 2)(2iV + L + 1) ^{2L - 5)(2L - 4)(2L - 2)(2L)(2L + 2)(2L + 3) 
(2iV- l)(2iV) (2L + 1)(2L- l) 2 



(2JV + L + 1) y/^iV — L)(2N + L + 3) (2L - 1)(2L)(2L + 4) 
(2JV-1)(2JV) (2L + l)(2L + 3) 2 
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(2N - L)y/(2N -L + 2)(2N + L + 1) (2L - 2){2L + 2)(2L + 3) 
(2JV-1)(2JV) (2L + 1)(2L-1) 2 



The only approximation made in the derivation of this equation in Appendix 3 has been the 
replacement of the reduced matrix elements of the quadrupole operator by their values in the 
contraction limit, which contain the intrinsic quadrupole moment qo, according to Eq. (|60p . 
Again the terms d)s and s^d give equal contributions of the same sign to the final result, 
despite the fact that they appear in Eq. ( TTOj) with opposite signs. 

One can simplify Eqs. ( !T6|) and (Fl8|) by making the following approximations 

1) Ratios of terms containing 2N can be replaced by unity, since we work in the large N limit. 

2) Ratios of terms containing 2L can be replaced by unity, if L is not too small. We will show 
below numerically that this approximation is an accurate one. 

Using these approximations we obtain Ri = 1 and R2 = 16/3. Replacing then the intrinsic 
quadrupole moment go by its value (N\/2) from Eq. fl55l) . we see that the two matrix elements 
vanish in the large N limit if the condition 

X(V) = — V ~ — (19) 
X{V) 8 (I-77) 2 1 ' 

is fulfilled. 

We remark that this condition is very similar to Eq. ( fl3|) . since y/7/8 = 0.3307 appears m 
the former, while 1/(2^2) = 0.3536 appears in the latter. This can be visualized in Fig. 2, in 
which these two conditions are compared to the original expression for the arc, Eq. (JUJ). 

It should be noticed that Eq. ( fl9l) has been derived in the special case of considering 
matrix elements within the ground state band only. This is similar to the condition discussed 
earlier in the context of Eq. (|S]), where the Q^jr 3 -\ operators simplify when the ground state 
band becomes energetically isolated. 

Concerning the accuracy of the approximations made above, one can check numerically the 
ratio R2/R1, for which the value 16/3 has been used. For N = 250 (the boson number used 
in Ref. [9]) one can easily see that the exact values of R2/R1 deviate by less than 1% from 
the approximate value 16/3 for angular momenta between 8 and 88, while they deviate by less 
than 5% for angular momenta between 4 and 176. Therefore the approximation is accurate 
for a large fraction of the spectrum, the lower-lying one. In the present case of iV = 250 the 
ground state band extends up to L = 500, thus the approximation is good (deviations less 
than 5%) for the lowest 1/3 of the spectrum. 

It should be clarified that in both subsections 3.1 and 3.2 the intrinsic quadrupole moment 
is finally replaced by its value, given by Eq. (15 5|). The main difference, though, is that in 
subsection 3.1 the intrinsic quadrupole moment appears by approximating the quadrupole 
operator by the intrinsic quadrupole moment, by looking at the SU(3) Casimir operator, as 
explained in detail below Eq. (1551) . while in subsection 3.2 the intrinsic quadrupole moment 
arises through the proper detailed calculation of the matrix elements of the commutators of 
the relevant parts of the Hamiltonian with the quadrupole operator within the ground state 
band. Therefore subsection 3.2 serves as a detailed justification of the approximation made in 
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subsection 3.1 . This detailed justification has been carried out within the ground state band, 
which suffices in the present case, since we work in the contraction limit, in which the ground 
state band gets isolated from the rest of the bands. 

It should be further noted that the method used in the present subsection, i.e., the detailed 
calculation of matrix elements, leading to Eq. (1191) . allows for an estimation of the accuracy 
of the approximations used within the ground state band for each value of L, as discussed 
above for the N = 250 case. In contrast, the method used in subsection 3.1, leading to 
Eq. ([TBI , uses an approximation (the replacement of the quadrupole operator by the intrinsic 
quadrupole moment), which cannot be tested separately for each state. The same holds for the 
work of Alhassid and Whelan (2J [3] , leading to the regular region approximated [6 J by Eq. ( HUH) . 
Indeed, in Refs. [21 EI], the statistical properties of the spectrum as a whole are considered, for 
iV = 25. This implies that numerical studies for large boson numbers (N = 250, for example), 
using statistical tools, should be undertaken, considering the spectrum as a whole, as well as 
by parts (lowest 1/3, middle 1/3, highest 1/3, for example). Studies of this kind should further 
clarify the relation between the present results and the Alhassid-Whelan arc of regularity. 



4 The SU(3) symmetry 

The question is now raised about what happens in the triangle formed by U(5), 0(6), and 
SU(3) [lj, the algebra containing the quadrupole operators with x — +v7 /2, which is known 
to correspond to oblate nuclei, while SU(3) is related to prolate nuclei. 

It turns out that the relevant calculation follows the same steps, with two notable differ- 
ences: 

1) — yf/2 is replaced by +\/7/2 everywhere. 

2) The intrinsic quadrupole moment changes sign (see Appendix 2), in agreement to the 
well known fact that the intrinsic quadrupole moment has positive values for prolate nuclei 
and negative values for oblate nuclei [I]. 

As a result of these two changes, Eq. ( !T3|) takes in the large N limit the form 

= -W2 (T^j) + T- (20) 

It becomes then clear that for a given rj in this case x acquires the opposite value from the 
one it gets within the U(5)-0(6)-SU(3) triangle. This is in agreement with the well known 
fact that properties within the U(5)-0(6)-SU(3) triangle are mirror images of the properties 
appearing within the U(5)-0(6)-SU(3) triangle [TJ. 



5 The 0(6) symmetry 

The successful determination of a line in the symmetry triangle characterized by approximate 
SU(3) symmetry, raises the question if a similar line related to the 0(6) symmetry can be 
determined. The 0(6) algebra [TJ is generated by the angular momentum operators of Eq. 
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flU) and the operators (d^d)^\ forming together the 0(5) subalgebra, plus the quadrupole 
operators 

Qok^^d + Js)^. (21) 
Following the same procedure as in Section 3, we are going to examine the conditions un- 

" (2) 

der which the Hamiltonian commutes with the generators Q (e) £■ The needed commutation 
relations are listed in Appendix 1. 

The first term of the Hamiltonian gives 

[Hi, QgJ 6)i J = cn[n d , gg ( } 6)i j = cn[(dh - Jd)<?\ (22) 

which is similar to Eq. (]6]). 
Using 

Q { 2z = Qo m + x(Jd)f\ (23) 
in the second term of the Hamiltonian, one gets the intermediate result 

+x{0d)f\ Q^MU + Qi%dW-l Qol)J-}} (24) 
Again this expression can be simplified in the large N limit, since in this limit the contraction 
of 0(6) to R 5 [SO(5)] [22] [23] takes place (see Appendix 4). This means that the commutators 
in the first two terms will vanish. Since 

[Qom* Q%),»} = 2(202i/|l£ + ^d)fl + 2(2£2z,|3£ + v)(Jd)f} u , (25) 

the vanishing of this commutator implies that terms of the form (d^d)^ can be ignored. 
Eq. ([21]) can be rewritten, without using any approximations yet, as 



k=l,3 ' 

+ xii^s - Jd)^Qf)W + {Qf\d)s - S tJ)P))W]. (26) 

In the large N limit the terms containing (d)dp k > can be omitted. In addition, in this limit 
Qq\q\ can be replaced by the intrinsic quadrupole moment (a scalar), which is iV in the present 
case (see Appendix 4). Then in the large N limit one is left with 

[Qf ■ Qf, Q { oUJ = 2N X [(d ] s - Jd)f. (27) 

The commutator for the second part of the Hamiltonian in the large N limit then reads 

[H2, Q% ))V ] = -2c(l - r]) X [(dh - (28) 

In order to get a vanishing commutator, the coefficients of (d^s — s^d)^ in Eqs. ( 122]) and f ]28]) 
should cancel, leading in the large N limit to the condition 

XM = J^-y (29) 

However, this does not suffice yet to guarantee the existence of a line corresponding to the 0(6) 
symmetry. One has also to consider the commutators of the Hamiltonian with the generators 
of 0(5), to be considered in the next Section. 
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6 The 0(5) symmetry 



As it has already been mentioned, the 0(5) algebra is generated by the angular momentum 
operators of Eq. (0J and the operators (Sd)f\ We should examine the conditions under 

which the commutator of the Hamiltonian with (cftd)^ vanishes. 

In this case the first term of the Hamiltonian makes no contribution, since is known to 
be an 0(5) scalar [5]. Taking advantage of the fact that Q Q ^ • Qo{&) ls a ^ so an 0(6) scalar [5] 
in order to simplify the calculation, one obtains 

[H, (dtj)(3)] = ^ x [(( d tJ)(2) (s tJ +d t s )(2))(3) _ (( a tJ +d t s )(2) (d t^P))(3)] 
_^^10 x [(( d tJ)(4)( s tJ +d t s )(2))(3) _ (( S tJ +d t s )(2)(dtJ)(4))(3)] 

_ ^VlO X 2 [((dtJ)(4)( d tJ)(2))(3) _ ((dtj)(2)( d t d ~)(4))(3)]_ (30) 

In this case no simplification due to contractions can be made. In order to get a vanishing 
commutator one needs to put x = 0, thus being confined on the U(5)-0(6) side of the triangle. 
This finding is in agreement with the known existence of an 0(5) subalgebra along the U(5)- 
0(6) side of the triangle [5]. 

Going back to the question of the existence of a line related to the 0(6) symmetry, we 
see that in Eq. (1291) one has to put x = 0, as required by the 0(5) subalgebra. Then one 
ends up with rj = 0, which represents the 0(6) vertex alone. Thus no line related to the 0(6) 
symmetry exists within the triangle. 



7 A different parametrization 



In recent years, an IBA Hamiltonian in common use reads [241 EH 

c 



H((,x) 



AN, 



-Q { ? ■ Q? 



B 



(31) 



where c is a scaling factor, while the rest of the symbols have the same meaning as in Eq. 
([1]). The parameter ( ranges from to 1, while the parameter \ ranges from to -a/7/2, as 
above. The parametrizations of the Hamiltonians ([T]) and (131 p are related by [7] 



4(C~1) 
3C-4 



c=-a(4-3C). 



(32) 



From these relations and the results reported at the end of subsection 3.1 we find that for large 
iV the critical line at x — corresponds to ( crit = 0.5, while at x — —V7/2 it corresponds to 
C c „* = 8/17 = 0.471. 

Eq. (fl3|) . valid in the large N limit, in this parametrization reads 



X (C) = V2 --1 



V7 

2 ' 



(33) 
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while Eq. (fT9|) . again valid in the large N limit, reads 



X(0 = ^T 7-1 -^r = ^T 7-2 • (34) 



8 Discussion 



In the present work we find the first analytical evidence for the existence of an approximate 
symmetry inside the symmetry triangle of the IB A. The SU(3) symmetry found extends from 
the SU(3) vertex until close to the critical line separating the spherical and prolate deformed 
shapes/phases, following the Alhassid-Whelan arc of regular behavior amidst chaotic regions. 
Thus it also points to an underlying SU(3) symmetry as an explanation for the existence of this 
arc of regularity. The present line of SU(3) symmetry has been determined in the limit of large 
boson numbers, taking advantage of the contraction of SU(3) to [R 5 ]SO(3), the algebra of the 
rigid rotator, in this limit. The proof is valid for the lower part of the spectrum, since SU(3) 
irreps with A much higher than the angular momentum L have been used in the contraction 
procedure. The contraction of 0(6) to [R 5 ]SO(5), the algebra of the 7-unstable rotator, has 
also been worked out, leading to a similar line, which however is shrinked to a point when the 
0(5) symmetry is imposed. 

The approximate symmetry determined here bears some similarity to the concept of quasi- 
dynamical symmetries (QDS) [TTI [T2l [13j [T31 [15], in which several features of a symmetry 
persist far outside its expected region of applicability. A difference is that in the present ap- 
proach an analytical proof of the existence of the symmetry is offered, while QDS have been 
so far located numerically. 

A concept related to the present work is that of partial dynamical symmetries (PDS) 
[29| 1301 |3T| [32] , situations in which part of the states preserve all the dynamical symmetry 
(PDS type I), or all the states preserve part of the dynamical symmetry (PDS type II), or part 
of the states preserve part of the dynamical symmetry (PDS type III). The present approach 
can be seen as an analog of a PDS of type I, since part of the states (the low lying ones) 
preserve all the dynamical symmetry. 

The method used in the present work is of wider applicability. It can be used for alge- 
braic Hamiltonians known to possess symmetries due to existing subalgebras, materialized 
for specific parameter values. The present method allows the extension of the region of the 
parameter space of the Hamiltonian in which the specific symmetry appears, in cases in which 
the relevant subalgebra contracts to a different algebra in some limiting case (the large boson 
number limit in the present case). 

The application of the present techniques to IBA-2 [1], in which distinction between protons 
and neutrons is made, to the sdg-IBA model [TJ [33], in which the L = 4 boson is taken into 
account, to the spdf-IBA model [311, ESI [36], in which negative parity bosons with L — 1, 3 
are included, and to the schematic Hamiltonian of Ref. [221 EZJ EB], which exhibits a PDS, are 
interesting tasks. Shape/phase transitions have recently been studied in both the parameter 
space of IBA-2, which is a tetrahedron [391 HQl E], and the parameter space of sdg-IBA [12], 
which is a prism. 
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Appendix 1. Commutation relations. 

We list here the commutation relations needed for obtaining the results of Sections 3-6. They 



are obtained through standard angular momentum coupling techniques |45j. 

[n d , (efts + Jd)f] = (Ss - Jd)f\ (35) 

[{Sd)f\ (Ss + S td)W] = (2£2z/|2£ + v)(Ss - s^d)fl u , (36) 

[h d ,(rfd)¥ ) ]=Q, k = 0, 1,2,3,4, (37) 

\{Sd)f ) {Sd)^] = -lQ^2i2vU + u)\ 2 2 I l}(Jd)fl, (38) 

[( d t s + jtjjW {d \ s + s tJ)(2)] = 2 £ ( 2^|^ + v)(dtd)?l, (39) 

fc=l,3 

[(j 8 + s tj)f , (dtj)«] = -y|( 2 av|2e + i/)(<ft a + s^lfi, (40) 

[(Ss + S tJ)f , (dtj)(3)] = _^| (2 £3z,|2£ + + s^lfi, (41) 

[(Sd)f\ (rft^W] = __L(i^|i e + ^(rftj)^, (42) 

[(^ ) ,(^«] = 2Vl5(2eii/|2e + i/){2 2 2} (dtJ) S- (43) 

[(dtj)f,(e^ 2 2}^+"- (44) 



Appendix 2. The SU(3)^ [R 5 ]SO(3) contraction 

The SU(3)— > [R 5 ]SO(3) contraction has been studied in Refs. [191 BO] - It is a procedure 
in which the full SU(3) algebra, consisting of 8 noncommuting generators, is shrinked into 
an SO (3) algebra (consisting of 3 noncommuting generators), accompanied by 5 mutually 
commuting operators (the quadrupole operators). This simplification occurs in the limit of 
large boson number in which, in SU(3), all intrinsic excitations rise in energy, isolating the 
ground state band so that SU(3) goes over, approximately, into a simple rigid rotator. The 
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resulting algebraic structure is, indeed, known |21j to be the algebra of the rigid rotator. The 
scaling of Ref. [1] is used here. (The quadrupole operator in Refs. [19j [20] is 2v2 times the 
quadrupole operator of Ref. pp.) 

The SU(3) commutation relations read 

[L 6 , L v ) = -V2(l£li/|lf + v)L i+v , (45) 
[4 Qsl( 3 ),v\ = -V&m^ + v)Q% 3U+u , (46) 



The second order Casimir operator is 

C 2 [SU(3)} = I 



?o (2) • n {2) + - f ■ f 

z VS[/( 3 ) Vsc/(3) 4 



(4f 



while its eigenvalues in the Elliott basis, (A,/x), are 

C 2 (A,/x) = ^(A 2 + /x 2 + A/x + 3A + 3/x). (49) 

If we consider SU(3) irreducible representations (irreps) with large values of C 2 (A,/x), that is 
for large boson numbers, we can rescale the quadrupole operator as 

Q {2) 

-(2) Vst/( 3 ),£ . . 

Qsu(3U = m A \ ( 5 °) 



The first two commutation relations remain unchanged by the rescaling, while the last one 
becomes 

= i/I (2e2z/|ie+z/) cfey (51) 

Then in the limit of large values of C 2 (A, /x) one gets 

[*lsu(3),z> ^su(3),J = 0- (52) 

This result, which is obtained for large boson number, is called the contraction of SU(3) to 
[R 5 ]SO(3), where [R 5 ]SO(3) is the algebra of the rigid rotator [21], generated by the angular 

(2) 

momentum operators of SO(3) and the five commuting operators g^™^, £ = —2, —1, 0, 1, 2. 

An immediate consequence of Eqs. ( l5Ti) and ( 1521) is that, in the contraction limit, terms 
proportional to the angular momentum L can be ignored. In the IBA framework, in which L 
is proportional to (Sd) 1 , as seen in Eq. (j3J), this implies that (d^d) 1 terms can be ignored. 

In the limit of large values of C 2 (A, /x) and A > /x the intrinsic quadrupole moments become 

MM 

i _ i 



g = ^=(2A + /x + 3), g 2 = - if) (/x + + 2), (53) 
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where K is the eigenvalue of the angular momentum projection on the body-fixed z-axis, for 
which K < L is valid, as one can see from the algorithm of the SU(3)dSO(3) reduction pQ. 
(Remember at this point that the quadrupole operator used in Refs. [HI EUJ [46] is 2\/2 times 
the quadrupole operator used in the present work.) For states with A >> L (therefore also 
A >> K) and A >> /i one then obtains [19] 



* = ^, (54) 

while q 2 becomes negligible. Since the ground state band belongs to the (2N, 0) irreducible 
representation (irrep) of SU(3), while other low-lying bands belong to irreps (2N — 4i — 6j, 2i), 
i = 0,1,2,. . . , j = 0,1,2,. . . with relatively low i, j, the contraction occurs in the large N limit. 
Thus in the case of interest the intrinsic quadrupole moment becomes 



qo 



NV2. (55) 



An equivalent statement is that one can approximately replace the operator Q s Ij^ by the 

scalar A/ a/2, as one can see from Eqs. ( 1481) and ( 1491) . since the terms containing L and \x 
are negligible in this limit, having as a consequence that only the first term in the rhs of 
these equations survives. A formal justification for this replacement is given in subsection 
3.2 and Appendix 3, where matrix elements of the commutators of the relevant parts of the 
Hamiltonian with the quadrupole operator are properly considered, resulting in the appearance 
of the intrinsic quadrupole moment. 

It should be noticed that the above results have been obtained in irreps with A >> L, thus 
they regard the low lying part of the spectrum. 

In SU(3) the irreps are built out of the (2,0) irrep, while in the case of SU(3) the irreps are 
built out of the (0,2) irrep pQ. As a result, in the SU(3) framework one is interested in states 
with large values of C 2 (A,/i) and A < fi, in which the intrinsic quadrupole moments become 

1 



qo = + 2/i + 3), q 2 = — A V3(A - K)(X + K + 2). (56) 

For states with fi » L and fi » A one then obtains 

qo = -± (57) 



while q 2 becomes negligible. Since the ground state band belongs to the (0, 2N) irrep of SU(3), 
while other low-lying bands belong to irreps (2i,2N — 4i — 6j), i = 0,1,2,. . . , j = 0,1,2,. . .with 
relatively low i, j, the contraction does occur in the large N limit, the intrinsic quadrupole 
moment becoming 

go = -NV2. (58) 



Since SU(3) is associated to prolate shapes, while SU(3) is related to oblate shapes, the signs 
in Eqs. (153]) and (1551) are consistent with the fact that intrinsic quadrupole moments are 
known to be positive for prolate nuclei and negative for oblate nuclei [I]. 
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Appendix 3. Matrix elements 



In order to show how Eq. (fTTjl is derived, we consider in detail the matrix element 
( [N] , (2N, 0) , x = 0, L 1 1 ( (rft a ) (2) Q g> (3) ) (2) 1 1 [N] , (2iV, 0), x = 0, L) 
= v / 5E{l 1 L 2 // }([iV],(2iV,0),x = 0,L||( C /t s )(2)||[ i v],(2iV,0) ) x = 0,L'') 

([iV], (2N, 0), x = 0, L"| |QS (3) | I [iV], (2JV, 0), x = 0, L), (59) 

where the standard formalism concerning the matrix elements of the tensor product of two 
tensor operators (13] has been used. By x we denote the Vergados quantum number [27] 
[not to be confused with the parameter x of the Hamiltonian (pQ)], which corresponds to 
an orthogonal basis, while the Elliott quantum number K [UJ, coinciding with the angular 
momentum projection on the body-fixed z-axis used in Appendix 2, corresponds to a non- 
orthogonal basis. 

This expression can be simplified since in the contraction limit the matrix elements of the 
quadrupole operator become jl6] 



([N], (2N,0),x = 0,L'\\Qf u{3) \\[N}, (2N,0), X = 0, L) = v / 2LTT(L020|L / 0)g , (60) 

where go is the intrinsic quadrupole moment of Eq. ( 153]) . This can be seen from the Elliott 
matrix element jUJ of the quadrupole operator. 

Using Eq. (ISP]) in Eq. ( 159]) one obtains in the contraction limit 

([An,(2A^0),x = 0,L||((d^ 
£ {I I } (^020|Z,"0)<[A^], (2A^, 0), x = 0, Z| | (rft s )(2) 1 1 [A r], ( 2Ar, 0), x = 0, L">. (61) 

The matrix elements appearing in the last equation can be written in terms of matrix elements 
of single boson operators 

({N},(2N,0),x = 0,L\\(dh)W\\{N},(2N,0),x = 0,L'') = V5{l l , \ 

([N], (2N, 0), x = 0, L\\d\\ [N-1], (2N - 2, 0), x = 0, L") 
([N-1], (2N - 2, 0), X = 0, L"\\s\\ [N], (2N, 0), * = 0, L"). (62) 

The matrix elements of single boson operators needed here and in subsequent calculations 
are listed in Table 1. Some of them are given in Ref. PQ, while the rest have been calculated 
following the method of Ref. jlH], using the triple-barred SU(3) reduced matrix elements given 
there and the SU(3)dSO(3) coefficients of Vergados [27] . 
The final result in the contraction limit reads 

([A^], (2iV, 0), x = 0, ^| |((^ S )^g^)^| | [iV], (2A^, 0), x = 0, i.) = 

l_ N pyr—r [2 (2N - L)(2N + L + 1) (2L - 3) 2 (2L + 5) 2 

14^2 3 (2N-1)(2N) (2L - 1) 2 (2L + 3) 2 
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(2N-L- 2)J(2N - L)(2N + L + 3) J(2L - 1)(2L)(2L + 2)(2L + 4)(2L + 6)(2L + 7) 



(2iV- l)(2N) 



(2L + l) 2 (2L + 3)< 



(2iV + L — l)y (2iV - L - 2)(2N + L + 1) ^(2^ - 5)(2L - 4)(2L - 2)(2L)(2L + 2) (21, + 3) 
(2N-1)(2N) (2L + 1)(2L-1) 2 

In the same way in the contraction limit one finds 

([iV],(2iV,0),x = 0,L||(g( 2 )(rft s )( 2 ))( 2 )||[iV],(2iV,0),x 



0,L) 



=N\/2L + 1 

14>/2 

(2iv + x + i; 



2 (2N - L)(2iV + L + 1) (2L - 3) 2 (2L + 5) 2 

3 (2N- 1)(2N) (2L- l) 2 (2L + 3) 2 

(2A - L)(2N + I + 3) (2L - 1)(2L)(2L + 4) 

(2N-1)(2N) (2L + l)(2L + 3) 2 

(2iV - L)MN - L + 2)(2iV + L + 1) (2L - 2)(2L + 2)(2L + 3) 



(2AT-l)(2iV) 



(2L + l)(2L-i; 



.(63) 



(64) 



For the terms involving s^d in the same way one finds 

([iV],(2iV,0),x = 0,L||(( S tJ)( 2 )g( 2 ))( 2 )||[iV] ) (2iV,0),x = ) L) = 
-([iV],(2iV,0),x = 0,L||(g( 2 )(rft s )( 2 ))( 2 )||[iV],(2iV,0),x = 0,L), (65) 
([iV],(2iV,0),x = 0,L||(g( 2 )( S td)( 2 ))( 2 )||[iV],(2iV,0),x = 0,L) = 
-([iV],(2iV,0) ) x = ) L||((rft s )( 2 )g( 2 ))( 2 )||[iV],(2iV,0),x = 0,L). (66) 

Using these results in the calculation of the matrix elements of [H2, Qsu(3)\ [including only 

1 we find Eq. 



the terms appearing in the second line of Eq. 



TUT), with replaced by Q 



SU(3)h 



Appendix 4. The 0(6)^ [R 5 ]SO(5) contraction 

A procedure similar to that of Appendix 2 is followed in the contraction of 0(6) to [R 5 ]SO(5) 
[22| [23] . This is a procedure in which the full 0(6) algebra, consisting of 15 noncommuting 
generators, is shrinked into an SO (5) algebra (consisting of 10 noncommuting generators), 
accompanied by 5 mutually commuting operators (the quadrupole operators). The resulting 
algebraic structure is known [23] to be the algebra of the 7-unstable rotator. 
The commutation relation for the quadrupole operators reads 

[Qolw Qoku\ = 2 E + ")(Jd)fl. (67) 

fe=l,3 

The second order Casimir operator is [1] 

C 2 [0(6)} = 2Q% } • gg ( } 6) + 4 £ (rftj)« . (dtj)« (68) 

fe=l,3 
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Its eigenvalues are 

C 2 (a) =2a(a + A), (69) 

where a is the quantum number characterizing the irreps of 0(6). 

If we consider 0(6) irreps with large a, we can rescale the quadrupole operator as 

42) _ Vo(6)U (7()) 
\JC 2 {(T) 

Then the commutation relation for the quadrupole operators becomes 

H&M. «o(e),J = 2 £ PC + -)^^#- (71) 

fc=l,3 L ^ a > 

Then in the limit of large a (and small t, where r is the quantum number characterizing the 
irreps of 0(5) ) [23] one gets 

[QoU,v 9o(6),J = °- (72) 
This procedure is called the contraction of 0(6) to [R 5 ]SO(5), where [R 5 ]SO(5) is the algebra of 
the 7-unstable rotator, generated by the operators of SO (5) and the five commuting operators 

<7o(6),f! £ = — 2, —1, 0, 1, 2, which are the coordinates [23] . 

An immediate consequence of Eqs. ( J7TI) and ( 1721) is that, in the contraction limit, terms 
proportional to ((Pd) k can be ignored. 

The most leading 0(6) irrep, to which the ground state band belongs, is (N). Thus in the 
large boson number limit it is appropriate to use this contraction. The intrinsic quadrupole 
moment will then be 

Qo = *, (73) 

as can be seen from Eqs. fl68l) and fl69|) . Thus in the case of interest the intrinsic quadrupole 
moment becomes 

Qo = N. (74) 

It should be noticed that the above results have been obtained in irreps with a » r, thus 
they regard the low lying part of the spectrum (since L < 2r, as seen from the algorithm of 
the SO(5)DSO(3) reduction [1]). 
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Table 1: Single boson matrix elements derived according to Ref. [35] . See Appendix 3 for 
further discussion. 



([N], (2N, 0), x = 0, L|| S t|| [iV - 1], (2iV - 2, 0), * = 0, L) = ^^^+^ ^Ny/2L + T 



([A], (2A, 0), x = 0, L||rft|| [A - 1], (2A - 2, 0), * = 0, L) = -^ j^^^^ ±|ViVV2T+T 
([A], (2 A, 0), x = 0, L||dt||[JV - 1], (2 A - 2, 0), * = 0, L + 2) = s^^^^^^^N^L + l 
([A], (2A, 0), x = 0, L| |dt| | [A — 1], (2A — 2,0),% — 0, L — 2) — ] /^^^^^VnV2LTT 



([A - 1], (2A - 2, 0), X = 0, L\\s\\ [A], (2A, 0), * = 0, L) = ~ A / ^gggjgff^ VW2L + 1 



([A - 1], (2A - 2, 0), x = 0, L|M | [A], (2N, 0), x = 0, L) = -/ g:g?ff^ ^v^ 



([A - 1], (2A - 2, 0), x = 0, L\\d\ | [A], (2A, 0), x = 0, L + 2) = f 2 ^^ v / iVv / 2T+5 
([A - 1], (2A - 2, 0), x = 0, L||d| | [A], (2A, 0), x = 0, L - 2) = ^ pE^^^i |VA V2L^3 
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Figure 1: (Color online) IB A symmetry triangle in the parametrization of Eq. ([1]) with the 
three dynamical symmetries, the Alhassid-Whelan arc of regularity [Eq. flHJ)], and the present 
line of Eq. ( 115]) (labelled as analytic). The shape coexistence region [TU] between spherical 
and deformed phases is shown by slanted lines near the U(5) vertex. In addition, the loci of 
the degeneracies E(2p)=E(2+) (dashed line on the right, corresponding to the SU(3) QDS 
discussed in Ref. [9]) and E{Ai)=E{^2 ) (dotted line on the left, also discussed in Ref. [9]) 
are shown for Nb=250 (top) and Nb = 25 (bottom). In the middle and bottom parts, the 
//-diagram and the A-diagram, based on Ref. [3j, are shown. See subsection 3.1 for further 
discussion. 
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Figure 2: (Color online) Location of the arc of regularity, as described by the original Eq. 
(fUj) . and as predicted by the findings of the present work, Eqs. ({TBI and (TT9l) . The r\ axis 
has been reversed, in order to correspond directly to Fig. 1. See subsections 3.1 and 3.2 for 
further discussion. 
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